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We study momentum imbalance as a function of jet asymmetry in high-energy heavy-ion collisions. 

To implement parton production during the collision, we include all Leading Order (LO) 2 —>■ 2 and 
2 —>■ 3 parton processes in pQCD. The produced partons lose energy within the quark gluon plasma 
and hadronize collinearly when they leave it. The energy and momentum deposited into the plasma 
is described using linear viscous hydrodynamics with a constant energy loss per unit length and a 
total energy loss given by a Gaussian probability centered around a mean value £ and a half-width 
AS. We argue that the shape of the asymmetry observed by the CERN-CMS Collaboration can 
indeed be attributed to parton energy loss in the medium and that a good description of data is 
achieved when one includes a slight enhancement coming from the contribution of 2 —>■ 3 parton 
processes that modifies the asymmetry distribution of the dijet events. We compare our results to 
CMS data for the most central collisions and study different values for S and AS. 

PACS numbers: 25.75.-q, 25.75.Gz, 12.38.Bx 


I. INTRODUCTION 

Strongly interacting hard probes are useful tools to 
study the properties of the quark-gluon plasma created 
in a relativistic heavy-ion collision. These hard probes, 
energetic quarks or gluons, are produced in the hard scat¬ 
tering that occurs at the earliest stage of the collision. 
They transfer energy to the medium via multiple scat¬ 
terings and eventually emerge as hadrons. The way the 
energy is transferred to the medium can be studied to 
infer some properties of the plasma such as temperature, 
viscosity, density, size, etc. To this end an important 
observable is the momentum imbalance between leading 
and away-side jets which allows one to investigate how 
momentum gets distributed in the away-side. Another 
such observable is the asymmetry distribution as a func¬ 
tion of the jet asymmetry Aj defined by 


PT,1 + PT,2 

where Pt,i i Pt ,2 the transverse momenta of the recon¬ 
structed leading jet and subleading jet respectively. The 
jet asymmetry is therefore positive definite by construc¬ 
tion. When both jets have similar momenta, Aj is close 
to 0 whereas a monojet (that usually corresponds to the 
case where the jet in the away-side is totally quenched) 
has asymmetry Aj = 1. 

The CMS collaboration has measured this asymmetry 
for jets with cone radius R = 0.5, triggering on leading 


jets with pT,i > 120 GeV and on subleading jets with 
Pt,i > 50 GeV [l|,0. Taking the case reported in Ref. [I|, 
for the most central PbPb events where the data set con¬ 
tains only jets with pT.i < 300 GeV, the asymmetry dis¬ 
tribution shows a significant deficit of events with Aj ~ 0 
as well as a significant excess of events with Aj ^ 0.4 as 
compared to a PYTHIA J- HYDJET simulation. This 
behaviour in the measured asymmetry distribution has 
been studied under different approaches: Monte Garlo 
simulations with improved jet finding algorithms, para¬ 
metric studies of energy loss models in connection with 
jet finding algorithms, models where NLO pQCD be¬ 
haviour of jet in-medium evolution is incorporated and 
solving relativistic hydrodynamic equations coupled with 
hadron momentum distribution profiles. The inter¬ 

play of jets and the hydrodynamic evolution of the bulk 
was first discussed in the context of Mach cone formation 
M and has been extended recently to other observables 
7|. The angular broadening of a medium-induced QCD 
cascade has also been recently studied in Ref. Q 

With these studies at hand, recent reviews on the state- 
of-the-art and future of jet physics in heavy ion collisions, 
attribute the measured jet asymmetry mostly to energy- 
loss mechanisms and argue that it is possible to go a 
step forward and look for other properties that could be 
constrained by this asymmetry (see for example [1,I3 and 
references therein). 

Another observable that can be measured as a func¬ 
tion of Aj is the missing transverse momentum. This is 




2 


constructed on an event by event basis, projecting the 
momentum of measured tracks onto the leading jet axis 
and summing over all tracks with transverse momentum 
larger than a minimum value p™™ >0.5. 

= X! ~Pt' (2) 


This quantity is negative definite in the leading jet 
hemisphere and positive definite in the away-side hemi¬ 
sphere 0,i. The missing momentum in the away-side 
can be obtained either for tracks within the reconstructed 
jet around the subleading hadron, in which case one ex¬ 
pects that the momentum outside the jet compensates 
the missing one inside the jet, or else, on the whole away- 
side hemisphere. The data in Ref. [J] shows that a large 
negative contribution to the average missing momentum 
(^1) (in the direction of the leading jet) for px > 8 
GeV, is balanced mostly by the combined contributions 
from the 0.5-8 GeV regions outside a cone of angular 
size A(p < 7r/6 centered along a direction opposite to the 
leading hadron. There is however also a small contribu¬ 
tion outside this cone from tracks with px > 8 GeV with 
0.3 <Aj< 0.5. 

The question we want to address in this work is 
whether these two observables, the asymmetry distribu¬ 
tion and the missing transverse momentum, can be better 
characterized if we account for the possible contributions 
from partonic processes where the hard collision of two 
partons in the incoming protons/nuclei gives rise to three 
partons in the final state. Qualitatively, two hard partons 
in the away-side that subsequently hadronize produce a 
broader asymmetry distribution and a harder component 
of the missing momentum outside the reference cone of 
the subleading jet which goes more in line with the re¬ 
ported measurements. This kind of studies can be im¬ 
portant to carry out a more precise determination of the 
energy loss mechanism. 

In two recent works [l3, [HI we have studied the effect 
of hadron production from 2 —>■ 3 parton processes in az¬ 
imuthal angular correlations using the leading order QCD 
matrix elements. These works have been extended to the 
description of the energy-momentum deposited into the 
medium [mil using linearized viscous hydrodynamics. 
In this work we apply this formalism to study the relative 
contribution from 2 —3 and 2 —>■ 2 parton initiated jets 
in the analysis of both the asymmetry distribution and 
momentum imbalance. The work is organized as follows: 
In Sec. m we collect the ingredients needed to describe 
the energy-momentum deposited in the medium by fast 
moving partons using linearized viscous hydrodynamics. 
In Sec, mil we present the results of this analysis and show 
that by including a 20% contribution from 2 —3 parton 
initiated processes one can achieve a better description of 
the asymmetry and momentum imbalance distributions. 
We finally summarize our results and conclude in Sec. lIVI 


II. SIMULATION 

Our procedure consists of generating either 2 —2 or 
2 —>■ 3 parton scattering events using Madgraph Cl. 
In a given event we let one of the partons hadronize 
in vacuum [ll which is then labeled as the leading 
hadron, whereas the other partons travel through the 
medium and lose energy before they hadronize outside 
the medium. 

Several energy loss scenarios have been considered to 
describe hadrons and/or jets emerging from a relativis¬ 
tic heavy ion collision. For example, in Refs. ii the 
authors make use of fragmentation functions, to help dis¬ 
tinguish the stage of the collision where a given energy 
loss mechanism dominates. In this work we avoid im¬ 
plementing a specific mechanism but instead consider a 
simple and general scenario that encompasses many fea¬ 
tures of the ones discussed in the literature. 

We consider a scenario where the hard scattering that 
gives rise to the final state partons happens near the 
medium’s surface, such that the leading parton goes 
into vacuum whereas the loosing energy partons travel 
through most of the medium. To allow for an uncer¬ 
tainty in the location of the hard scattering we suppose 
that AE is normally distributed according to a probabil¬ 
ity density function. 

We take the energy lost by a parton from a Gaus¬ 
sian probability density distribution p{£) with an av¬ 
erage £ = 50 GeV and a width AS = 10 GeV. This 
corresponds to a situation where the leading hadron is 
emitted near the surface. For definitiveness we consider 
that this happens at AL = 5/6 L, where L = 10 fm is the 
total medium’s length and that the energy-loss per unit 
length as dE jdx = 5 GeV/fm. A parton with energy E 
produced in a hard scattering will travel an average dis¬ 
tance AL in the medium and lose energy AE (which on 
average is f ^ L{dE/dx) = 50 GeV) to emerge from the 
medium with energy E — AE and fragment in vacuum. 
The above values emerge from a study of azimuthal par¬ 
ticle distributions implemented in the context of energy 
deposition into the medium using linearized viscous hy¬ 
drodynamics [l3 that produces particles by means of a 
Gooper-Frye formalism. We assume that the sub-leading, 
but otherwise fast partons, produced in the hard scatter¬ 
ing that travel through the medium do not change their 
direction of motion due to interactions with the medium. 
This seems to be a plausible assumption since the final 
hadrons, that these partons give rise to, are still very 
energetic. The interaction of the fast partons with the 
medium gives rise to energy and momentum deposited 
into the medium described by linear viscous hydrody¬ 
namics which in turn produces low momentum partons 
distributed around the fast ones according to the Cooper- 
Frye distribution. 

Therefore, to compute 5p, that is the average px de¬ 
posited by the away-side parton into the medium, we 


3 




FIG. 1. Simulated dijet asymmetry distribution for leading jets of pt,i > 100 GeV and subleading jets of pT ,2 > 30 GeV (left) 
and pt ,2 > 50 GeV (right). The histograms show pp 2 —> 2 (hne dotted line) and AA 2 —>■ 2 (dotted line) and 2—^3 (solid 
line), separately. 


write [13 

Ar(Aj/)2 

dpT—e^P[-PT/To\ 

Se gySin(/) + g^cos(/)'\ 

4eo eo(l + c2) J 

The integral over </> in Eq. m determines in- and out- 
of-jet cone regions where the average pt is distributed 
and Ay, At the rapidity and freeze-out time intervals, 
respectively. 

The energy density 6e and the components for the mo¬ 
mentum density gi = (gT)j + (gL)^ {i = y,z, T and L 
are the transverse and longitudinal modes) are given by 



,5e = 






Ise{a,P) 


(4) 


and 



^3[T/^ (^) 


where the integrals (3) are obtained from the 

solution of the linearized viscous hydrodynamic equa¬ 
tions with a source term corresponding to the current 
produced by a localized and fast moving parton which is 
in turn proportional to a constant ener gy l oss per unit 
length. These are given explicitly in Ref. [l^ and written 
in terms of the dimensionless variables 

a = \z- vt\/ and /3 = xp/ ’ (6) 

which represent the distance to the source jz —and the 
distance from the parton along the transverse direction 


XT (along the y axis in the geometry we are using) in 
units of the sound attenuation length 

r =_ ^ _ (7) 

*-3eo(l + c2)- 

As implemented in Ref. [l^ Cq is the static back¬ 
ground’s energy density and Cs the sound velocity. This 
background is perturbed by the energy deposited by the 
fast parton moving near the speed of light. We focus on 
events at central rapidity y ~ 0 and take the direction 
of motion of the fast parton moving in the away-side to 
be the z axis (with the beam along the x axis). With 
this geometry, the transverse plane is the y — z plane 
and therefore, the momentum four-vector for a (mass¬ 
less) particle is explicitly given by p^ = {E,px,Py,Pz) = 
(pt, 0,Pt sin (j),PT cos (p), where p is the angle that the 
momentum vector p makes with the z axis. 

Equation ([3]) can be used in a two ways: on the one 
hand, by integrating the azimuthal angle over the full 
away-side hemisphere one can equate the average mo¬ 
mentum deposited into the medium with the total en¬ 
ergy lost by the fast parton and therefore extract the 
characteristic time At for the process. In this way we get 
At ~ 8.5 fm. On the other hand, given a certain amount 
of energy lost extracted randomly from the gaussian pro¬ 
file, by restricting the azimuthal angular integration, one 
can compute the amount of momentum going into the 
in- and out-of-cone regions around the away-side leading 
hadron. 

The CMS collaboration supplemented their study with 
an analysis of the momentum imbalance inside and out¬ 
side a cone with a hxed radius, as a function of the jet 
asymmetry Aj Qi- The analysis was done accounting 
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FIG. 2. Simulated dijet asymmetry distribution for leading jets of pt,i > 100 GeV and subleading jets of pt ,2 > 30 GeV (left) 
and pt ,2 > 50 GeV (right) The histograms in the top row show the benchmark asymmetry distribution for the pp 2 —>■ 2 sample 
compared to GMS data for pp collisions in the 0-10% centrality bin. The histograms on the bottom row show the asymmetry 
distribution for the AA 2 —>■ 2 (solid line) 2 —>■ 3 (fine dotted line) and combined (dotted line) samples, compared to GMS data 
for Pb-Pb collisions in the 0-10% centrality bin. The combined sample has 80% of 2 —>■ 2 and 20% of 2 —>■ 3 AA events. 


for the different px contributions in each Aj bin. The re¬ 
sults showed that most of the momentum in the leading 
side is balanced with momenta in the 0 to 8 GeV/c range 
in the away-side, collected within a cone whose angular 
distance from the leading hadron was A^i _2 > 27r/3. 

To account for the contribution of soft particles pro¬ 
duced by the hadronization of the momentum deposited 
into the medium by the fast parton Sp in our model, we 
compute the in-medium jet asymmetry variable as 

Aj = (8) 

Pt.i + PT,2 

where pt ,2 = Pt .2 + dp, when we are balancing the mo¬ 


menta inside the cone, with dp given by Eq. (|3]) and 
Pt ,2 = AE — dp, when we do it outside the cone. In the 
next section we look at observables plotted as functions 
of Aj as defined in Eq. © for pp, or as functions of Aj 
as defined in Eq. m for AA. Since what matters is the 
momentum of the final hadron and not the momentum 
of the initial parton used to compute the asymmetry, we 
simply call the asymmetry variable Aj. 
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FIG. 3. (Color online) Average missing transverse momentum for leading jets of pt,i > 100 GeV and subleading jets of 
Pt ,2 > 50 GeV, as a function of the event asymmetry. The left (right) column shows the balance for pp (AA) events arising 
from 2 —>■ 2 (top row) and 2 —>■ 3 (bottom row) partonic events. All panels show the negative momentum imbalance in the 
direction of the leading jet is shown above the horizontal axis and the momentum deposited into the medium is below the 
horizontal axis in the away-side hemisphere Acf) > n/2. 


III. RESULTS 

In this section we present the results of the asymmetry 
and momentum imbalance analysis on the simulated sam¬ 
ples, and compare this analysis to CMS data as reported 
in references [J, . We focus first on the asymmetry dis¬ 

tributions for the 2—^2 and 2 —>■ 3 proton-proton and 
nucleus-nucleus samples together with a combined sam¬ 
ple. The combined sample has 80% of 2 —>■ 2 prong events 
and 20% of 2 —>■ 3 prong events for both pp and AA. 

Figure [1] shows the simulated dijet asymmetry distri¬ 
bution for leading jets of px.i > 100 GeV and subleading 


jets of pt .2 > 30 GeV (left) and pt .2 > 50 GeV (right). 
The histograms show the benchmark asymmetry distri¬ 
bution for the pp 2 —> 2 sample (fine dotted line), the AA 
2 —2 sample (dotted line) and the AA 2 —>■ 3 sample 
(solid line). Note that the pp 2 —>■ 2 simulation is made 
out of symmetric events with Aj in the range 0 to 0.3. 
On the other hand the AA 2 —J- 2 data shows a depletion 
of symmetric events and an enhancement of asymmet¬ 
ric ones. The shift in the peak and the widening of the 
distribution from pp to AA data, is purely due to the en¬ 
ergy loss of the away-side parton that was implemented 
to generate this sample. The figure also shows a mutu¬ 
ally complementary effect in the AA sample between the 
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FIG. 4. (Color online) In-cone average missing transverse momenta for leading jets of px.i > 100 GeV and subleading jets of 
Pt,2 > 50 GeV, as a function of asymmetry. The figure shows the in-cone balance for AA events arising from 2 —>■ 2 (left panel) 
and 2 —>■ 3 (right panel) partonic events. The negative momentum imbalance in the direction of the leading jet is shown above 
the horizontal axis and the momentum deposited into the medium is below the horizontal axis, inside a cone of A<j) > 37r/4 
(7? = 0.5). 


2 —> 2 and the 2 —3 results: in the most symmetric bins 
{Aj < 0.3) the 2^2 events contribute more whereas in 
the most asymmetric ones {Aj > 0.3) the 2 —> 3 events 
give a larger contribution as expected. The histograms 
on the left panel have a wider distribution than those 
on the right one which also have more symmetric events. 
This is caused by the lower cuts on px on the left side 
panels compared to more restrictive ones for the panels 
on the right. 

Figure [2] shows simulated asymmetry distributions for 
leading jets of px.i > 100 GeV and subleading jets of 
Pt ,2 > 30 GeV (left) and pt ,2 > 50 GeV (right). In the 
upper row we plot the asymmetry distribution for the 
benchmark pp 2 —>• 2 sample and compare to CMS data 
for pp collisions in the 0-10% centrality bin. The lower 
row shows the asymmetry distribution for the AA 2 —2 
sample (solid line), the 2^3 sample (fine dotted line) 
and the AA combined sample (dotted line) compared to 
CMS data for Pb-Pb collisions in the 0-10% centrality 
bin 0,i- There is a small difference in the description 
of the data when using the combined AA sample rather 
than the A A 2 —> 2, on the central asymmetry bins. The 

2 —>■ 3 sample generates a slight depletion of symmetric 
events and a similarly small excess of asymmetric events. 
The interplay between this sample and the 2 —>■ 2 can 
account for some of the behaviour attributed to pure e- 
loss. The fact that we have a sample with 20% of events 
in a 2 —>■ 3 configuration and are still able to describe the 
data supports the idea that a dijet sample with true 2 —>■ 

3 initiated dijet events can be a component of the analysis 
that should be understood before coming to conclusions 
about the energy-loss mechanism at play. 


The asymmetry distribution analysis can be comple¬ 
mented with a study of the average missing momentum 
as a function of jet asymmetry. This can help visualize 
the effect of the energy deposited into the medium for 
different cone radii. In our analysis it is crucial to see 
how much of the missing momentum for a given cone ra¬ 
dius can be a attributed to the third hadron in a 2 to 3 
event. Using Eq. m we can construct histograms show¬ 
ing the effect of energy loss both in the whole away-side 
hemisphere and only within a cone of size tt/2 centered 
around the leading hadron, in the away-side. 

Figure [3] shows the average missing transverse momen¬ 
tum for leading jets of pt,i > 100 GeV and subleading 
jets of pt ,2 > 50 GeV, as a function of the event asym¬ 
metry. The left (right) column shows the imbalance for 
pp (AA) events arising from 2 —>• 2 (top row) and 2 —>■ 3 
(bottom row) partonic events. All panels show the nega¬ 
tive momentum imbalance in the direction of the leading 
jet (red histograms) and the momentum deposited into 
the medium (blue histograms) in the away-side hemi¬ 
sphere A(j) > 7r/2. There are two main features worth 
noting in these plots: First, the pp sample fills up to 
Aj = 0.5 but the AA sample fills up to Aj = 0.7. There¬ 
fore there is a magnification on (p!|,) when going from pp 
to AA (both in the 2 —>■ 2 and 2 —J- 3 samples) and this is 
due exclusively to the energy loss mechanism. Second, for 
the AA sample, the higher the Aj, the more important 
the 2 —>■ 3 component becomes with respect to the 2 —>• 2 
sample, with up to 30% effect on the most asymmetric 
bins. 

We now study the amount of missing px contained in- 
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FIG. 5. (Color online) In-cone average missing transverse momenta for leading jets of px.i > 100 GeV and subleading jets of 
Pt,2 > 50 GeV, as a function of asymmetry. The panels show the in-cone balance for AA events arising from 2 —>■ 2 (shaded) 
with (pJ|) scaled to 80% and the 2 —>■ 3 (solid) sample with (p!]l) scaled to 20%. The sample used was obtained with the energy 
lost by a parton taken from a Gaussian distribution probability density p{£) with an average f = 30 GeV and a half-width 
AS = 3 GeV (a) and with S = 50 GeV and half-width AS = 5 GeV (b). The negative momentum imbalance in the direction 
of the leading jet is shown above the horizontal axis and the momentum deposited into the medium is below the horizontal 
axis, inside a cone of At)) > Svr/d {R — 0.5). 


side a cone of a given radius, as a complement of the half¬ 
hemisphere analysis done in Fig. |31 Figure 0] shows the 
in-cone average missing transverse momentum for leading 
jets of pt,i > 100 GeV and subleading jets of pt ,2 > 50 
GeV, as a function of asymmetry for AA events arising 
from 2 —>■ 2 (a) and 2 —s- 3 (b) partonic events. Shown 
are the negative momentum imbalance in the direction 
of the leading jet and the momentum deposited into the 
medium inside a cone of > 37r/4 {R = 0.5). Note 
that the relative contribution of the in-cone momentum 
on the away side is smaller for the 2 —?► 3 as compared to 
the 2 —>• 2 case. This can be understood as a consequence 
of the second particle in the away side not being emitted 
within the reference cone. We emphasize that the analy¬ 
sis shown in Figs. [3] and 21 agrees with the recent study 
of jet formation and bulk evolution interplay of Ref. Q. 

The GMS collaboration has reported an analysis 
whereby the balance of momentum is achieved by pro¬ 
jecting all tracks onto the leading jet axis. From this 
sample a further selection of in- and out-of-cone tracks is 
made to study how much momentum is pushed outside 
the away-side cone, presumably by energy loss effects. 
We now make a comparison with GMS data and learn 
about the role the 2^3 processes have in the energy 
loss mechanism, which in turn will give us a rough es¬ 
timate of their contribution to (px). The comparison 
can only be done with balanced events at the hadronic 
level, so we implement a fragmentation where momentum 
is almost conserved (the fraction of hadron momentum 


from the parent parton is almost the same in the lead¬ 
ing and the away sides), to emulate the fact that GMS 
has already projected all tracks onto the leading jet axis. 
Figure O shows the average missing transverse momen¬ 
tum for leading jets of pt,i > 100 GeV and subleading 
jets of pt ,2 > 50 GeV, as a function of asymmetry, aris¬ 
ing from this AA sample, compared to GMS data [I|. 
The histograms correspond to the 2^2 (shaded) sam¬ 
ple with (pl|.) scaled to 80% and the 2 —>• 3 (solid) sample 

with (p!],) scaled to 20%. The AA sample used was ob¬ 
tained with the energy lost by a parton taken from a 
Gaussian probability density distribution p{£) with an 
average f = 30 GeV and a half-width AS = 3 GeV for 
panel (a) and with f = 50 GeV and half-width A£ = 5 
GeV for panel (b). Note that a smaller average energy 
loss and a narrower width are more consistent with the 
data. Also, note that the most asymmetric bins in both 
panels are exclusively populated by the 2 —>■ 3 sample, 
since the 2^2 sample cannot reach this asymmetry 
bin. This is due to the fact that with a balanced sample 
the asymmetry in the 2 —s- 2 case comes exclusively from 
energy loss. Therefore for the values £ ± AS considered 
and the implemented px cuts, the largest asymmetry bins 
cannot be populated. This is not the case for the 2^3 
sample, since one of the away side partons, could end up 
outside the reference cone. 







































IV. SUMMARY AND CONCLUSIONS 

In this work, we have studied the way that fast moving 
partons deposit their energy and momentum when travel¬ 
ling in a medium. We perform this study by looking into 
the momentum imbalance as a function of the jet asym¬ 
metry, using generated 2 —>• 2 and 2^-3 parton events 
that lose energy and hadronize collinearly to form jets. 
To estimate the average energy loss deposited into the 
medium by the fast moving partons, we used linearized 
viscous hydrodynamics and the Cooper-Frye formula to 
calculate how this energy is distributed in different mo¬ 
mentum contributions in the away-side. We argue that 
for conditions resembling those achieved in heavy-ion col¬ 
lisions, the shape of the obtained asymmetry and momen¬ 
tum imbalance agree with the ones reported by the CMS 
Collaboration. This in turn can be explained as origi¬ 
nating from the energy loss of the partons created in the 
hard scatterings and travelling in the plasma, including 
a slight enhancement produced by the contribution from 
2 —)• 3 events. 

Furthermore, we found that the (p!j<) observable shows 
that the contribution of 2 —3 events is enhanced when 
going from proton-proton collisions to nucleus-nucleus 
collisions, with up to 30% effect on the most asymmetric 
bins and agrees with the ones presented in recent studies 


of jet-bulk interplay Q. Finally, we perform an analysis 
with a balanced hadron momentum sample of the aver¬ 
age missing pT to compare with CMS data. We show 
that this analysis favours a smaller average energy loss 
and a narrower width, compared to the one used for our 
asymmetry distribution analysis. 

All together, our results suggest that an analysis con¬ 
taining a mixed two- and three-jet sample with a realistic 
energy-loss profile may prove useful for a better charac¬ 
terization of the jet momentum imbalance as a function 
of jet asymmetry. Progress in this direction will be re¬ 
ported elsewhere. 
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